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ABSTRACT 

We derive analytically the vorticity generated downstream of a two-dimensional rippled hydro- 
magnetic shock neglecting fluid viscosity and resistivity. The growth of the turbulent component of 
the downstream magnetic field is driven by the vortical eddies motion. We determine an analytic 
time-evolution of the magnetic field amplification at shocks, so far described only numerically, until 
saturation occurs due to seed-field reaction to field lines whirling. The explicit expression of the ampli- 
fication growth rate and of the non-linear field back-reaction in terms of the parameters of shock and 
interstellar density fluctuations is derived from MHD jump conditions at rippled shocks. A magnetic 
field saturation up to the order of milligauss and a short-time variability in the X-ray observations of 
supernova remnants can be obtained by using reasonable parameters for the interstellar turbulence. 
Subject headings: Physical Data and Processes: turbulence; ISM: cosmic rays, magnetic fields 
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1. INTRODUCTION 

Compelling evidence has been cumulated that individ- 
ual shell-type Supernova Remnant (SNR) shocks accel- 
erate charged particles, i.e., electrons and probably ions, 
up to energies at le ast of the order of 10 13 — 10 14 eV (see 
e.g. Cassiopeia A (jAharonian et al. || 2001[ ), R X J1713.7- 
39 46 jAharoirian et al.||2004| ), Tycho's SNR ( |Acciari et 
aL"||20fT )). Charged particles are likely to be accelerated 
by two simultaneous mechanisms: the so-called Fermi 
first-order, i.e., repeated shock crossing of the particle 



( Axford et al.||1977| |B ell 197 8a|b| [Bland Ibrd «V Oinkn 

ire at 



978| iK rymskii 1977j), and drift along the shock surface 
I Jokipii 1982, 1987). From detection of non-thermal X 



ray rims ( |Vmk and Laming 2003 Bamba et al.| [20041), 
rapid time- scale variability of A -ray hot spots dUc hiyama 
et al. | [2007l) and 7-ray emission in extended regions ( |Ac-| 
ciari et al.|2011 ), a magnetic field at the shock far exceed- 
ing the theoretically predicted shock-compressed field 
has been inferred. The X- ray rims could be due to damp- 



ing of the magnetic field rtPohl et al. 2005 ) rather than 



to synchrotron emission, although corresponding narrow 
filaments in the r adio emission have not been observed 
( jRothe nflug 20 04]) ( see ho wever high-resolution radio im- 
ages in[Dyer et al. (2009)). Magnetic field amplification 



might be also relevant to in situ measurements of the 
plasma downstream of the solar- wind termination shock 
( Burlaga et al.|2007[ ), where fluctuations have been mea- 
sured of the same order as the mean, or to radio observa- 
tions of Mpc scale sh ocks at the edge of galaxy clusters 
( Briiggen et al.|[201~2 ). Whether or not such a magnetic 
field amplification m SNR is to be associated with ener- 
getic particles at the shock is still subject of controversy. 
So far, SNR observations could not rule out ci- 
ther of the two following mechanisms of field ampli- 
fication: a) microscopic plasma instabilities generated 
by cosmic-rays current flowing upstream of the s hock 
and therein exciting non-resonant magnetic modes (Bell 
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20041; b) macroscopic turbulent fluid motion down- 
stream of the shock seeded by inhomogeneities of the 
upstream medium triggering vortical eddies and tangling 
the magneti c field lines hence am plifying the turbulent 
component (Giacalone & Jokipii 2007) (cosmic-rays pres- 
sure gradient has also been propose d as driver of the am- 
plification in iDrury & Falle (1986)). The former mech- 
anism was also discussed in the kinetic theory approach 
( Amato fc Blas i 2009); however, nume rical simulations 
(see, e.g., Riquelme & Spitkovsky ( 2009 1 ) could find only 
a moderate amplification (B/Bq ~ C(10)), and the un- 
folding of its non-linear extension and observational im- 
plications is currently active. 

Collisionless shocks propagate in turbulent and inho- 
mogeneous media undergoing rapid corrugation of their 
ideal planar surface. Numerical simulations have shown 
that the picture of a planar shock is inappropriate to 
describe secular evolution of downstream medium. The 
unshocked medium might be strongly inhomogeneous at 
several scales and the cold interstellar clumps strongly 
deform the shock surface. The passage of an oblique 
non-relativistic shock through inhomogeneous medium 
has been known fo r longtime to generate vorticity in the 



downstream flow (Ishizuka et al. 1964); in a conduct- 
ing fluid the turbulent motion at scale I with fluid ve- 
locity vi and local density p leads t o an exponentially 



ampl ified magnetic field B 2 = Airpvf (Landau & Lifshitz 
19601. Such a dynamo action amounts to a systematic 



conversion of the fluid kinetic energy into magnetic en- 
ergy at each scale separately, possibly until equipartition 
is reached ( |Kulsrud||2005[ ). Recent numerical 2D-MHD 
simulations have shown that such an amplifi cation can 



be very efficient (Giacalone & Jokipii 
numerical studies using the thermal ins' 



2007). Further 
lability of the in- 



terstellar medium (ISM) turbulence, i.e., condensation 
of t he inter stellar gas due to catastrophic radiative cool- 
ing (Field 196 5J), confirme d the efficient magnetic field 
growth (( jlnoue et al.|2012 ) and references ther ein). Two- 
di mensio nal simulations of relativistic shocks ( Mizuno et 
aT[|2011 ) show that small-scale dynamo can operate also 



downstream of the shocks of Gamma-Ray Bursts out- 
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flows, suggesting that the dynamo action downstream of 
shocks might shed light on the energy equipartition at 
magnetized shocks. 

In this paper we provide an analytic derivation of 
the vorticity generated by clumpy unshocked medium, 
hence of the magnetic field amplification, downstream 
of a non-relativistic rippled collisionless shock. We ap- 
ply the Rankinc-Hugoniot jump conditions locally down- 
stream of an MHD shock to compute the vorticity gen- 
erated downstream. For the sake of simplicity, a two- 
dimensional shock is considered, i.e., observables depend 
only upon two space coordinates. The downstream vor- 
ticity depends on the magnitude of the tangential com- 
ponent of the velocity (shear), although a different in- 
terpretation is given here in terms of density gradient 
at clumps boudary, and the curvature of the shock sur- 
face as pre viously foun d for pure ly hydrodynamic shocks 
( |Truesdell||l952|JKevlahan||1997[ ). We also compute the 
back-reaction to vortical motion of the seed magnetic 
field advected downstream, so far accounted for only nu- 



meric ally (e.g. [Giacalone fc Jokipii 2007 Inoue et al. 
2012|. Using the small-scale dynamo theory we deter 



mine the time-evolution of the turbulent magnetic field 
in the downstream fluid until saturation epoch. 

The encounter of a shock surface with a density clump, 
also called Richtmyer-Meshkov (RM) instability ( |Brouil-| 
lette 2002 1 , has been extensi vely investigated in plasma 
labora tory experiments (see |Dimonte fc Ra maprabhu 
(2010) and references therein^ JNumerical simulations 
of magnetic shocks proved that RM instability drives 
transie nt events in several re gions of the Earth magneto- 
sphere ( Wu fc Roberts|1999|). Recen t plasma laboratory 
experiments ( |Kuramitsu et al.|2011 1 made use of laser to 
test the magnetic held amplification by densi ty inhomo 



geneities at shocks of supernova remnant (see Sano et al. 
(2012) and references therein) 



1'he small-scale fluid vortices close behind the shock 
grow on time-scale smaller than the particle accelera- 
tion time-scale, which depends on the seed magnetic field 
orientation, the isotropy of the turbulent component of 
magnetic field and the dependence of the spatial diffusion 
coefficients, parallel and perpendicular, on the particle 
energy. Therefore, the vortical field growth is unaffected 
by the presence of cosmic-rays at the shock. Magnetic 
field may also be enhanced by fiel d line stretching due to 



Rayleigh- Taylor (RT) instability (|Jun et al.|1995[ ) at the 
interface between the ejecta and the interstellar medium, 
i.e., far downstream of the shock. In contrast with the 
vortical turbulence, late-time RT turbulence might be 
affected by the highest energy particl e gyrating in the 



down stream fluid far from the shock (Fraschetti et al 
20101. However, RT structures are unlik ely to reach out 



the blast wave (( |Fraschetti et aL 2010) and references 
therein) and therefore to interact with vortical turbu- 
lence. Thus the dynamo amplification local behind the 
shock can be temporally and spatially disentangled from 
the field line stretching due to RT instability. 

This paper is organized as follows: in Sect. [2] we intro- 
duce the constitutive equations and describe tne features 
of a rippled shock. In Sect. [3] we define the vorticity 
in the local rotated frame, compute the vorticity down- 
stream of the shock and interpret the result in terms of 
vorticity growth and field back-reaction. In Sect. [4] we 
apply the small-scale dynamo theory to determine the 



time-evolution of the turbulent magnetic field. In Sect 
[5] we identify the dependence of field growth and non- 
linear field back-reaction on the physical shock parame- 
ters and discuss the implications for recent X and 7-ra; 
observations of non-relativistic SNR shocks. In Sect 
we summarize our findings. 

2. RIPPLED SHOCK 

We consider the propagation of a 2D non-relativistic 
shock front in an inhomogeneous medium. The time evo- 
lution of the two driving independent observables in ideal 
MHD approximation, i.e., the fluid velocity v = v(x,y) 
and the magnetic field B = B(x,y), is given, with no 
viscosity or heat conduction and for infinitely conduc- 
tive fluid, by 

<9 t v + (v • V)v + — + -L[Bx (VxB)] = (1) 



p Airp 

P = Vx(vx B) 



(2) 



where p, P, are respectively density and hydrodynamic 
pressure of the fluid (here dt = d/dt) . Equation (T2| does 
not incl ude any field-ge nerating term, such as Biermann 
battery ( Kulsrud]|2005[ ) , as the fluid on both sides of the 
shock is embedded in a pre-existing magnetic field. Note 
that the current density carried by cosmic-rays is here 
neglected: we aim to identify the growth of the mag- 
netic energy as generated by the vortical motion of the 
background fluid only. Thermal dissipation reduces the 
energy deposited in the magnetic turbulence and will be 
considered in a forthcoming publication. 

The shock is a dynamic surface due to the interaction 
with the upstream clumps. The kinematics of a 2D hy- 
drodyna mic rippled shoc k propagating in inhomogeneous 
medium ( jRavindran fc Prasad| [l993) comprises of a se- 
quence of unstable configurations of the shock surface. 
The conservation of energy requires that sections of the 
corrugated shock surface being ahead or lagging behind 
readjust to smooth-out growing corrugations. The time- 
evolution of the shock surface is acco unted for by the 
rate change along the moving surface (Prasad 2001) of 
i?(i, x,y), i.e., the local angle between the average direc- 
tion of the shock motion and the local normal to the 
shock surface (see Fig[l]). We assume here that such 
a self-deformation process of the shock surface occurs 
on a time-scale much greater than the turn-over time of 
the smallest eddies in the downstream flow of the shock. 
Thus the shock profile is "frozen" during the exponen- 
tially fast amplification which proceeds as the field is 
advected downstream. 

3. VORTICITY DOWNSTREAM OF MHD SHOCK 

In the inviscid approximation used here the conserva- 
tion of the vorticity flux applies: the vorticity shock- 
generated is transported along the flow "frozen" into the 
fluid, as a consequence of Hclmholtz-Kelvin theorem. 
The medium upstream of the shock has zero vorticity. 
The vorticity is calculated downstream at a distance from 
the shock large enough that the shock is infinitely thin, 
i.e., the thickness of the shock is much smaller than the 
local curvature radius at every point of the shock surface. 

At a rippled shock the MHD Rankine-Hugoniot jump 
conditions cannot be applied globally as the directions 
normal and tangential vary along the shock surface. In 
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Fig. 1. — Encounter of a shock surface with density enhance- 
ment regions: forward and lagging behind regions are formed that 
generate vorticity in the downstream fluid. 

a 2D shock propagating at average in the direction x 
(Fig IT]), from the velocity field of the flow v = (v x ,v y , 0), 
the vorticity is given by u> = V x v = (0, 0, u) z ) and the 
component in the direction z, outgoing from the paper, 
by uj z — d x v y — d y v x (all quantities are independent on 
z and we used d Xi = d/d Xi ). We use a local natural co- 
ordinate system (n, s), where n — (cosz?(i, s),sim9(£, s)) 
is the coordinate along the normal to the shock surface, 
s = (sim9(£, s), — cosz?(i 1 s)) is the coordinate parallel to 
the shock surface (FigjlJ). Note that $ = $(t, s) varies 
only along the shock surface and not along the orthogo- 
nal direction n, thus is independent on n (a lso, the time- 
evolution of $ in Ravindran & Prasad (19931, Eqs. [2.21]- 
[2.23], does not depend on n). In the local frame (n, s) 



the z-component of the vorticity becomes 
lu z = d s v n - d n v s + v s d s tf, 



(3) 
,cos$ 



where v n — w^cosf? + Wj,sim? and v s — t^sim? 
are respectively the local component of the velocity flow 
normal and tangential to the shock (see Fig IT]). The local 
derivatives here can be expressed in terms of derivatives 
in directions (x, if) as d n = n ■ V = cos?9 d x + sim? d y and 
d s = s ■ V = sini? d x — cos$ d y . The vorticity generated 
downstream of the shock in the presence of a magnetic 
field B = (B n ,B s ,0), for lj z = upstream (detailed 
computation is in the Appendix), results in 



5oj, 



1 



1 



pC r l r 
d s SB 2 B n SB 



[C r } u d s (pC r ) 



8- 



4tt 



cW 



-in 



5\d s B n - d n B s 



(4) 



where C r = C — v n is shock speed in the fluid frame (C is 
the shock speed in the normal direction n), r = Pd/Pu is 
the compression ratio at the shock and Sf = [f]d — [/]« 
indicates the jump of / from upstream (u) to downstream 
(d). 

In the present paper, we deal with the hypothesis that 
the turbulent magnetic component is only generated in 
the downstream flow through dynamo mechanism, i.e., 
the turbulence is not required on both sides of the shock 
as in the diffusive shock acceleration. The contribution 
to strong dynamo amplification in the downstream flow 



from far-upstream fluctuations of magnetic and flow ve- 
locity fields on several scales will be included in a forth- 
coming work. If the magnetic field is quasi-perpendicular 
to the average direction of shock motion, charged par- 
ticles can be still efficiently accelerated to high-energy 
by the motional electric field through drifting along the 
shock, until they are advected in the downstream flow or 
escape upstream, provided they are not scattered back 
to the shock to further acceleration. In an oblique or 
quasi-parallel configuration, if the field is amplified only 
in the downstream flow as described here, other up- 
stream turbulence or pre-existing magnetic instabilities 
are needed to scatter energetic particles back across the 
shock, increase the residence time in the upstream region 
and release larger energy accelerated particles. Since no 
upstream magnetic turbulence is included here, at the 
present stage this model reconciles dynamo amplification 
with particle acceleration occurring in the same region 
of the shock only for quasi-perpendicular magnetic field. 
Turbulence of upstream medium allows an extension to 
oblique and quasi-parallel cases. 

We consider a seed-magnetic field upstream uniform 
and normal to the average direction of motion (B = 
(0,£>n,0), or B n = Bosiwd and B s = —Bqcosi), see 
FiglTb. For the first term in the second line of Eq. Q, 

we can write d s SB 2 = d s 5B 2 s 2B 2 5(sm2$cLi)) = 0, as 

(5(sini?) = 8(d s -d) = 0. The last term in Eq.Q vanishes: 
<J[B co8i?(a a i?)] = and d n B s = 0. 

Assuming that amplification is efficient at the smallest 
scales (see Sect. EJ, B n and SB S are drowned out by 
the turbulent components: SB S ~ —tBq — Bl urb + Bq ~ 
_gturb an( ^ £^ ^ gturb f or a perpendicular field (see 
also last paragraph in the present Section) . Therefore the 
factors SB S and B n in Eq.Q include both the impulsive 
shock compression and the turbulence amplification. We 
can conclude that the vorticity produced downstream of 
a 2D shock propagating in an inhomogeneous medium 
with a uniform perpendicular upstream magnetic field 
(same as for parallel shock as shown later) can be recast, 
neglecting obliqueness, in a simple form: 



\Su z \ 



1 



d s p + d s C T 



B n 5B s 

AnpCr 



d s $, (5) 



where SB S is the jump across the shock of the magnetic 
field in the direction locally tangential to the shock sur- 
face including the Rankine-Hugoniot compressed seed 
field and the turbulently amplified field and B n is the 
component in the direction locally normal to the shock 
surface including the unchanged Rankine-Hugoniot and 
the turbulent component. 

The factor (r — l)/r shows the dependence of Suj z on 
the shock compression. The first term in Eq.pJ contains 
the density gradient (d s p), hence having the role of the 
baroclinic generation of vorticity, i.e., the term VP x Vp 
appearing in the equation for u> time-evolution: when 
an upstream clump (clump size is assumed greater than 
shock thickness) crosses the shock the clump is locally 
heated with a maximum compression at the front of the 
cloud and smaller along the side of the cloud, whereas 
the densi ty gradient is dir ected across the fluid-clump 
interface ( Klein et al.p 994): thus, the gradients of den- 
sity and pressure are no longer parallel and vorticity is 
generated in the transition layer (Vp 7^ 0). The ther- 
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mal instabil i ty model for ISM predicts a broad range 
( |Field||1965| |Begelman fe McKee|[T990| for the transi- 
tion layer between the cloud and the ISM shocked gas 
as a function of the thermal conduction(see Sect. p)h. 
The second, or corrugation, term in Eq.([5]) results from 
the finite curvature radius of the shock, i.e., for a planar 
shock d s C r — at every point of the shock surface. As 
shown in the following Section, the purely hydrodynamic 
terms (baroclinic and corrugation) drive the small-scale 
magnetic field growth. 

The back-reaction of the small-scale turbulent field is 
represented by the last term: vortical eddies shear and 
whirl field lines around enhancing the turbulent field as 
long as the entailed magnetic tension grows to strength 
large enough to halt such a growth. Similar dependence 
of the field back-reaction on the turbulent-field Alfven 
speed can also be reasonably derived by dimensional ar- 
guments ( jKulsrud|2005[ ) but it follows here from the ap- 
plication of local Kankine-Hugoniot jump conditions at 
rippled shocks. Note that in the absence of shock rip- 
ples, i.e., for infinitely large local curvature radius, the 
term d s d — > and back-reaction of the amplified field 
on the vorticity becomes negligible. Also, the vorticity 
downstream of a shock wave in Eq.([5h holds regardless 
the speed of the shock wave (Alfven, fast or slow magne- 
tosonic speed). We recall that in contrast with 3D turbu- 
lence, the shear of the fluid velocity along the vorticity, 
i.e., (a; • V)v, vanishes in a 2D flow and the vorticity 
generation in a 2D fluid might be under-estimated. 

Note that if the magnetic field upstream is uniform 
and parallel to the average direction of motion (Bo = 
(Bq , 0, 0), or B n = B cost? and B s = B smd), the mag- 
netic term in Eq.([5| is unchanged. Therefore, assuming 
efficient turbulent neld amplification, a difference in the 
growth of |w|, i.e., in the saturated turbulent field, is not 
expected between the perpendicular and parallel field up- 
stream cases, c onfirming previous numerical findings (In- 
oue et al.|2012 ) (see also Sect. Eh. The reason is that the 
turbulent field isotropically becomes much greater than 
Bo, regardless the orientation of the seed-field upstream. 

4. TURBULENT FIELD AMPLIFICATION 



bulcnt field obeys an unbounded exponential amplifica- 



As a result of dynamo action (Kulsrud 2005), the 
growth of a turbulent magnetic field is known to be gov- 
erned by the fluid vorticity at each scale. The unper- 
turbed field is initially too weak to affect the fluid veloc- 
ity field and the turbulent field grows exponentially fast, 
i.e., on the time scal e of the smallest eddies turnover time 
( Balsara et al.|2005 ), until the magnetic energy produces 
non-negligi ble effects on the velocity field and then sat- 
urates. 

The time-evolution of the total magnetic field (EqJ2h , 
i.e., both ordered and turbulent components, is deter- 
mined here by using the small-scale dynamo theory ( |Kul-| 
srud|20 05). The vortical swirling exponentially amplifies 
also the mean field. However, since the amplification 
time-scale is of the order of the smallest eddies turnover 
time, the saturation occurs much faster at small-scale. 
No statistical assumption is made on the magnetic field, 
in contrast with the usual numerical approach. This is a 
relevant simplification as the spatial diffusion of charged 
particles in a turbul ence depends on the magne tic power 
spectrum (see, e.g., Fraschetti & Jokipii (2011 1). 

The small-scale dynamo theory predicts that the tur- 



tion at a rate j3 (Kulsrud 2005 
1992J): de/dt = 2/3e, where e 



|Kulsrud fc Anderson 
B 2 /8irp is the total mag- 
netic energy per unit of mass, including seed and tur- 
bulent fields. As shown in Kulsrud (2005), the isotropy 
and homogeneity of the velocity correlation entails the 
following simple relation between the amplification rate 
of e and the vorticity generated downstream of the shock: 
f3 ~ (7r/3)<5u; z . However, such a version of the dynamo 
theory cannot resolve the dilemma posed by the un- 
limited magnetic growt h and several non-linear versions 
have been elaborated ( Kulsrud||2005 ). 

We present in this paper the back-reaction of the field 
to the vorticity growth and the saturation of B at realis- 
tic rippled shocks. We argue that Eq.([5]) for \Suj z \ allows 
to explore the regime of self-controllea growth wherein 
the turbulent field becomes so strong to affect the veloc- 
ity field. The time-evolution of the turbulence amplifica- 
tion until the saturation epoch can be determined: if we 
recast Eq.([5]) as |<5w z | = (?>/it)(t~ 1 — ae), then e satisfies 



de 
df 



= 2(t~ 1 - ae)e 



(6) 



o»G r 



the local 



f r -r [{C r /p)ud s p + a s (j r \ is 

growth rate of e and a — (2ir/i)d s , d/C r is the lo- 
cal back-reaction; the initial condition for Eq. (ph is 
e(0) = e = v\/2 = B$/8irp. In Eq.Q we have 
assumed that the turbulence dominates over B , i.e., 
SB s /y/8irp ~ y/s and B n /y/8irp ~ y/e: the turbulence 
grows isotropically downstream at the shock curvature 
scale as a consequence o f the isotropy of the flow veloc- 
ity field ( |Kulsrud||2005] ). We note that the ISM clump 
size is much larger than the shock thickness and the shock 
crossing is not impulsive on the vorticity generation time- 
scale; thus, e in the bracket of Eq.fph is continuously gen- 
erated during the shock crossing and is not constant in 
time. 

Neglecting the time dependence of r (the magnetic 
mode s grow slowly for initially weak field (Kulsrud 
20051), the solution is readily found: 



£o 



■(*) 



B 



(t) 



2t/r 



\-aT(\~e 2t l^)v 2 A l2' 



(7) 



for a uniform average interstellar matter density. Equa- 
tion ^j\ provides the first analytic time-evolution of a 
self-controlled amplification of magnetic energy through 
fluid vortical motion generated downstream of a 2D rip- 
pled shock. 

For t <^i t the term proportional to v A in the denom- 
inator is negligible and B grows exponentially (B/Bo ~ 
e*' r ). For times <>t the term proportional to v\ dom- 
inates over the first term and the saturation is attained: 



B/B ~ ^2/{arv\), 



(8) 



corresponding to an increase B/B ~ Ma, where Ma = 
C r y/4:Trp/Bo is the seed- field Alfven Mach number (see 
Sect. ph. Note that this result applies only if the as- 
sumption that the turbulence B turh downstream domi- 
nates the compressed seed field tBq. For interplanetary 
shocks (Ma < 10) the measured downstream turbulence 
is rarely amplified to values much greater than the mean 
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Fig. 2. — Saturation of the total magnetic field for various shock 
speed C r is shown: CV = 1,500 km/s (M A = 50), CV = 5,000 
km/s (M A ~ 170), C r = 15,000 km/s (M A = 500), assuming 
_R C = 10 17 cm and l F = 10 16 cm, that results in r < £j?/C r ~ 3 
years for CV = 5, 000 km/s (i? = 0.1 rad, r = 4 and u. 4 = 10~ 4 c). 
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Fig. 3. — Saturation of the total magnetic field for various R c is 
shown, for C' r = 4, 800 km/s (M A ~ 160), and using « F =5x 10 16 
cm (i9 = 0.1 rad, r = 4 and Uj4 = 10 — 4 c). 



field ( Burlaga et al.|2007 1 , and the term rB is no longer 
negligible with respect to B turb and a different analysis 
is needed. 

5. DISCUSSION 

Whether the magnetic field is amplified by waves ex- 
citation in the shock precursor or downstream through 
the coupling of the eddies motion with the magnetic field 
lines, or a synergy of the two mechanisms, or other pro- 
cesses, is not observationally settled yet. Only the former 
has been largely explored in recent years. This paper 
presents for the first time an analytic approach to the 
latter mechanism. 

In what follows we estimate r and a showing that 
Eq.([7]) applies to the time-evolution of magnetic energy 
in young and middle-aged SNRs. We do not aim here at 
predicting a typical length-scale of the X-ray rims due 
to the large variety of SNRs taxonomy, but rather we 
show that reasonable values of the clump/fluid transi- 
tion layer, shock curvature radius and shock speed are 
capable to produce the inferred efficient magnetic field 
amplification. 
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Time (year) 

Fig. 4. — Saturation of the total magnetic field for various values 
of l F for C r = 4, 800 km/s (M A = 160) and R c = 3 x 10 17 cm 
(1? = 0.1 rad, r = 4 and v A = 10~ 4 c). 

5.1. Field growth 

From the discussion below Eq.([5]), the magnitude of 
the vorticity downstream \Su z \ is expected to depend on 
the thickness of the layer of the interface clump/fluid 
and not on the size of the clump. Thermal instability 
model of two-phase fluid suggests that such a thickness 
is given by the Field length £p (|Field||1965| |Inoue et 



(Field 
al. 2006). However, £ F dep ends on the 



reat conduction 
"( |Begelman fc McKee||1990[ ) and is thus uncertain as the 
thermal conductivity might be sensitive to the magnetic 
field. The previous history of formation and evolution of 
the supernova progenitor might also affect the properties 
of the surrounding medium, i.e., uniform or stellar wind 
density profile. 

The growth time-scale can be recast as r _1 = 
% T T L [(Cr/p)ud.p + d.C r ] ~ ^CriR- 1 + r F l ). Nei- 
ther of the terms -R"" 1 , £ F l can be neglected because no 
particular assumption has been made on R c , expected 
indeed to be comparable to the size of the clump that 
corrugated the shock, with respect to £p. From the pre- 
vious argument it follows: 



r 1 RJ F 
r-lCrR c + £ F ' 



(9) 



The growth rate r _1 of the turbulent field, and of 8u> z , 
increases with shock speed and it depends only on hy- 
drodynamic quantities, except a possible unknown de- 
pendence on B in £p here disregarded. If £p <IC R c , from 
Eq.@ it holds r - £ F /C r (see also FigJ4j>. Thus, the 
fieldamplification saturates faster in a region of smaller 
If providing a constraint testable by multiwavelength 
observations. This justifies the description in Sect|3]that 
the vorticity, and the strong turbulent field, are gener- 
ated in the transition layer separating the cold clump 
from warm ISM and does not depend on the clump size. 
Note that r is smaller than the travel-time of the 
shock across the interior of the ISM clump, i.e., cloud- 
crushing time t cc , resulting fro m the ram-pressur e equi- 
librium at the clump boundary ( Klein et al.|1994 ): t cc — 

\/Pc/poL/C r , where p c and po are respectively mass- 



G 
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density in the clump and in the surrounding fluid. Struc- 
tures of scale L stable under external heating or radiative 
cooling condense if L 3> £p; thus t cc — \J ' p c / 'p L/C r ^> 

y/pc/po^F/C r ~ ^p c /poT > t. Thus, for any p c /p the 
ISM clump survives across the shock provided that the 
density gradient at the clump boundary is steep enough. 

5.2. Non-linear back-reaction 

As the magnetic field strengthens, it reacts to field 
lines whirling halting the turbulence growth. In more 
general terms, as the field increases by dynamo effect it 
also releases its tension by unwinding at a rate of or- 
der of Alfven speed: the ba ckreaction grow s with the 
turbulent field Alfven speed ( Kulsrud 2005 1 . The local 
back-reaction of the field a ~ (J s t)/C r can be estimated 
by 

a ~ $/{R c C r ). (10) 

Clearly a is enhanced by small curvature radius (a ~ 
1/Rc) as the more corrugated is the shock, the shorter 
the eddies turnover time, the more efficiently the turbu- 
lent field is created downstream and back-reacts to fluid 
motion. From Eq.Q, it is clear that as e increases, the 
back-reaction term grows faster than the driving term, 
eventually becoming dominant: the strongest field has 
the smallest growth (cfr. the backreaction dependence 
on B in the non-linear Landau damping). 

5.3. Secular evolution of the turbulent field 

Fig{2] depicts the growth of the turbulent field for var- 
ious shock speeds, assumed constant in timer] given an 
ISM field of the order of Bo ~ 3/iG, the turbulent field 
saturates at B ~ 1.2-3. mG for C r = I, 500-5, 000 km/s 
on the year time-scale. Such a rapid growth of magnetic 
energy is compatible with X-ray obser vations of SNRs 
RXJ 1713.7 - 3946 (q < 4,500 km/s jUchiyama et jjT 
2007|) and Cas A (|Patnaude fc Fesen||2"0T)9[) brightness 



Ma- From Sect. [4] the saturation can be written as 



variations detected on year time-scale in small-scale hot 
spots structures, attributed to synchrotron electron cool- 
ing. Using R c = I0 17 cm and lp — I0 16 cm, we find an 
amplification to B ~ 3. mG within 3 years. Such a value 
of ip is to be compared with the spatial scale of the Chan- 
dra RX J17I3.7— 3946 bright spots, estimated as < 0, 03 
pc. Similar length (~ 1 14 — 1 16 cm) and time (~ 1 yr) 
scales are found in simulations of the effects of magnetic 
field turbulence on the obse rved synchrotron e mission 
images and spectra in SNRs ( |Bykov et al.||2008 |. Thus, 
the magnetic energy increase and the X -ray variability 
might have a time-scale (1 yr) much lower than the SNR 
hydrodynamic time-scale and might occur in middle-aged 
SNRs, not necessarily youn g SNRs (RX J1713.7 - 3946 
age is estimated as 1, 600 yr ( Stephenson & Green|: 
The high shock speed C r ~ 15,000 km/s in 
comparable to observations of the youn gest SNR in our 
galax y, i.e., 100 years old G1.9 + 0.3 (Reynolds et al. 
20081. Thus, a rapid field saturation even up to B ~ 10 
inGis predicted at SNR shocks within a few months. 

The computation above indicates a linear increase of 
the saturation value with the shock speed (see FigJ2J) , or 

2 On the year time-scale the speed C r is approximately constant: 
in the ejecta-dominated phase (t < 1, 000 yr) of a self-similar SNR 




The equipartition between the magnetic pressure of the 
turbulent component downstream of the shock, the ther- 
mal pressure and the fluid ram pressure (B 2 /8tt ~ pC 2 ) 
implies: (B/B ) 2 ~ (87r/Bg)(pCj?) ~ 2M A . Thus the 
scaling B/Bq ~ Ma indicates that the amplified field 
pressure might locally become comparable with the ram 
pressure (and the gas pressure), according to the factor 
(R c -\- lp) / (lp{f): the non-linear back-reaction, found us- 
ing jump conditions at ri ppled shock s in this paper, is 



no longer negligible (cfr. Bell (2004)) and the rates of 



growth and unwinding are equa 

Figl3] shows the turbulent field growth for various val- 
ues oiR c , with a shock speed C r = 4, 800 km/s (close to 
the highest measured for the SNR1006 shock (Katsuda 
eTaI]| 20l2))) : Fro m HST observations of SN1006 ([Ray 
a 



mond et a. 



20071, R c is inferred to be less thanT7TU 



of the forward shock radius, estimated as R s ~ 9.1 pc 
= 2.8 x 10 19 cm. However, the overall smoothness of 
SN1006 suggests that the ripples scale might be lower 
(R c ~ 10 16 cm) in more clumpy-structured SNRs. This 
justifies the smaller value of R c used in Fig(2] Moreover, 
the absence in SN 1006 of hot spots with rapid time- 
variability, as detect ed in synchrotron coo ling regions of 
RX J1713.7 - 3946 ( |Uchiyama et al |2007| , might result 
from the spatial location at high Galactic lat itude where 
sma ll-scal e clumps are not present (see also Katsuda et 



aLj | |20T0l ). Thus, a saturation B/B ~ 500 - 800 for 
R c = 10 iY c m is compatible with the optical constraints 
on SN1006 (Raymond et al. 2007), indicating that a re- 



gion with strong amplified field variable on 10 yrs time- 
scale could still be observed in future high-spatial resolu- 
tion efforts before the advected amplified field is drowned 
by projection effects. 

Despite several uncertainties, a numerical estimated 
range for lp in the I SM is given by [IP 1 6 — 10 18 1 cm 
(seeEq.s [5.1], [2.7] in |Begelman fc McKee| ( |1990| ». The 
field growth is clearly more efficient at the small scale 
If , as also shown in Fig|4j In FigJ3] an increase in R c 
reduces the field back-reaction (a ~ l/i? c )resulting in a 
larger saturation value. In contrast (Fig|4| an increase 
in £ F (Up <C R c ), reduces the growth rate (r ~ £ F /C r ) 
and therefore B/Bq. 

ISM density turbulence has been known for long-time 
through radio scintillation to obey the Kolmogorov scal- 
ing oyer several decades in len gth-scale (Lee and Jokipii 



1976 Armstrong et al. 1995) 



mj ect ion scale, L c 



from the outer scale, or 
If the 



contact discontinuity with radius R ~ t 
changes by a factor < 5 — 6. 



the speed CV ~ t 



BTTpc ( |Noutsos|[2012 | 
turbulent fluid velocity downstream of the shock follows 
a ID Kolmogorov power spectrum at a cert ain scale fc -1 
as nu merical simulations seem to suggest (Inoue et al. 
20121 (P v (k) ~ fc -5 / 3 ), the vorticity spectrum at scale 
k is given by P u (k) — k 2 P v (k) ~ fc 1 / 3 : in the inertial 
range the vortical energy at small-scale is greater than 
at large scale. Also, since the field at small-scale grows 
and saturates faster, at that scale t he magnetic en ergy 
coincides with the turbulent energy ( Kulsrud]|2005 ). We 
indicate here such a scale to be < 0.01 pc. However, 



as a result of the uncertainties on the Field length ip 
( Field| 19*65 ), within the scenario presented here a unique 
length-scale of magnetic growth cannot be identified. A 
range of scales consistent with the ISM fluctuations is 
shown here to account for optical and X-ray observa- 
tions. The confirmation will be sought through kinetic 
analysis or 3D numerical simulations in forthcoming pub- 
lications. 
An attractive feature of the microinstabilities mecha- 



nism in Bell ( 2004 1 is that magnetic field growth is asso- 
ciated with energetic particles at the shock and the field 
saturation is driven by the cosmic-rays of highest energy, 
possibly 10 15 eV protons. In contrast, the lack of direct 
association between particle acceleration and field ampli- 
fication in the dynamo mechanism here does not disfavor 
it as an interpretation of the X-ray rims in SNRs: the 
particle acceleration might proceed through the motional 
electric field E = |v x B|/c along the rippled shock sur- 
face (such an accelerat ion is known to be faster than in 
the parallel field case (Jokipii 1987)). We note that the 
scale of the smallest eddies, or highest turbulent field, in- 
ferred here (£p ~ 10 16 cm) is comparable to the gyroscale 
r g of energetic protons with energy E = 10 15 eV in a field 
B = 300^iG (r g = 10 16 cm); hence, the mean free path A 
of the highest energy particles, except the Bohm scatter- 



ing case, i.e.. 



A ~ r g , is typically greater than lp. If this 
happens, the turbulence downstream is able to enhance 
scattering of highest energy particles, thereby leading to 
an additional particle acceleration at the shock possibly 
beyond the knee of the cosmic-ray spectrum. 

6. CONCLUSION 

We have investigated the generation of vorticity down- 
stream of a non-relativistic two-dimensional rippled 
shock front typical of shell-like supernova remnant ex- 
panding in a turbulent interstellar medium. By using 



only jump condition at a rippled shock surface, we have 
derived the temporal evolution and the saturation of the 
turbulent magnetic field downstream of the shock, in- 
cluding the non-linear field back-reaction. We conclude 
that the saturation of B by small-scale dynamo action de- 
pends on the shock speed, on the thickness of the layer 
of the ISM clumps and on the shock curvature radius. 
The saturation value is found counter-intuitively not to 
depend on the size of the IS M clumps, as also found for 
purely hydrodynamic shocks ( Klein et al.|1994 ). Perpen- 
dicular and parallel field upstream cases lead to the same 
saturation values, if amplification is very strong, due to 
fast isotropization of the downstream turbulence. 

Our finding shows that small-scale dynamo might ex- 
plain non-thermal X-ray and optical observations of 
young/middle-aged supernova remnant shocks. The sec- 
ular evolution of the turbulent magnetic field derived 
here might help to shed light on the evolution of young 
SNRs to be discovered. The youngest SNRs in our galaxy 
or other "historical" SNRs might be in the growing phase 
of magnetic field described here. Current and next gen- 
eration of hard X-rays observatories (NuSTAR, Astro-H) 
will provide a helpful probe for the mechanism here de- 
scribed. 
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APPENDIX 



A. COMPUTATION OF THE VORTICITY 



The u) z generated downstream of the shock is found, for convenience, by computing separately the jump across the 
shock of every term on the right side in Eq.([3|, with the shock profile being frozen. Recalling that if B ^ the 
transverse flow velocity is not conserved acrossthe shock (5v s ^ 0), we find: 



Su z = d s (Sv n ) - S(d n v s ) + (8v a )d a <&. 



(Al) 



We compute t he term S(d n v s ) by using the equations of motion, i.e. Eqs.|l][2]), and following the algorithm used in 
Kevlahan (1997). The magnetic field B = (B ni B s ,0) ^ accounts for the back-reaction to the vorticity growth (see 
Fq.([5])). Rewriting the two components of Eq.pJ) in coordinates (n, s), one finds two new equations: 



costf T - 
sini? T - 



- sim? G = 
cos$ Q = 



(A2) 



where we have defined 



dv n d-d d n P B s 

F= —rr + v s — + h - — (d n B s - d s B n - B s d s -d) 

dt dt p 4np 



g= dv ± _ v d-d | d s P B n 
dt dt p 4:Ttp 



(d n B s -d a B n -B a d a ti), 



and d/dt = dt + v n d n + v s d s 



Fraschctti 



Equation ( A2 ) implies T = Q 



0. We recast Q 

did 1 dv s 

dt v n dt 

B n 

Anpv n 



= as (see also Sect. 2 in Ravindran & Prasad (19931) 

d s P 

pv n 

{-d n B s + d s B n + B s d 8 $) = 0. 



(A3) 



We now evaluate the jump of Eq.( A3) across the shock (8(d s d) = 0). Multiplying Eq.(A3l by (pv n )d and making use 
of the mass conservation (8[p(v n — C)] — — 5[pC r ] — 0), we obtain 



(pv n )dS(d n v s ) + 5(pd t v s + pv s d s v s ) ~ (pv n ) d (5v s )d s & 
+ d s SP + C8p(d n v a ) u - C6p(d t #) u - C6p(v a ) u d,0 



+ ^S[~d n B s 

47T 



d s B n + B s d s $] = . 



(A4) 



Note tha t if B = this expression is equivalent to Eq.s [2.8-2.10] in Kevlahan (1997) with d n d = 0. 

In Eq.(A4) the mass conservation readily gives: (pv n )d(5v s )d s -d + Cdp(Vs) u ds'd = 8(pv n v s )d s '&. We make here 
the assumption that the shock velocity in the tangential direction is much smaller than in the normal direction, 
i.e., obliqueness effects are neglected; hence the jump in transverse momentum writes S[p(v n — C)v s — B n B s /Air] = 
0, implying S(pv n v s ) = C8{pv s ) + B n 8B s j 'An. Differentiating with respect to s the jump condition of the normal 
momentum, i.e., 8[p{C r ) 2 + P — B 2 /8n] = 0, provides 



- d s SP = -d s [pC r Sv n ] + d s 6[B 2 /8ir}. 



(A5) 



We replace the expression for d s 5P from Eq .(|A5[) into Eq.(|A4|) and substitute the resulting expression (|A4|) for 
S(d n v s ) into the vorticity downstream (Eq. (lAll) 



Suj z 



1 

pC r 

d s 5B 2 

8tt 



By using the jump condition 5v n = - — -[C r ]«) one can write 
r-1, 



PV S 



Dv 

lot 



+ 2 BnSBs d s # 

An 



B n 

An 



'[C r ] u d s (pC r ) 



S[d s B n - d n B s 



(A6) 



where (Dv/Dt) = (dt + v s d s + C5 n )v. No t e that , if B = 0, the second line in Eq.(A6) vanishes and the remaining 
terms co rres pond to the result in |Kevlahan ( J1997 ). 

In Eq.(A6) we can approximate to first order in xf the shear term: v s - (Dv/Dt) ~ C r (D-d / Dt). The mass conservation 
S(pC r ) = for the first term in Eq.(A6) gives S[pv s ■ (Dv/Dt)] ~ pCrSlDtf/Dt] ~ pC r <5(u s )9 s i?. The time-derivative 
dd /dt ~ Dd/Dt for an M HD shock wa ve depends on shock speed, magnetic fi eld c omponents and density gradient 



in a complicated manner (Prasad 20011. Remarkably, in the shear term of Eq.( |A6[ ) the only non- vanishing jump is 
(Sv s )d s 'd thus independent on the particular equation of motion for i9. 

U sing the MHD Rankine-Hugoniot conditions for the transverse velocity, Sv s — —B n 5B s /(AirpC r ), we can recast 
Eq.(A6) as Eq.Q in the main text. 
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